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Abstract — We present a novel method for the evaluation of 
shot noise in quasi one-dimensional field-effect transistors, such 
as those based on carbon nanotubes and silicon nanowires. The 
method is derived by using a statistical approach within the 
second quantization formalism and allows to include both the 
effects of Pauli exclusion and Coulomb repulsion among charge 
carriers. In this way it extends Landauer-Biittiker approach by 
explicitly including the effect of Coulomb repulsion on noise. We 
implement the method through the self-consistent solution of the 
3D Poisson and transport equations within the non-equilibrium 
Green's function framework and a Monte Carlo procedure for 
populating injected electron states. We show that the combined 
effect of Pauli and Coulomb interactions reduces shot noise in 
strong inversion down to 23% of the full shot noise for a gate 
overdrive of 0.4 V, and that neglecting the effect of Coulomb 
repulsion would lead to an overestimation of noise up to 180%. 

Keywords - Shot noise, FETs, nanowire transistors, carbon 
nanotube transistors. 

I. Introduction 

In the last few years, a huge collective effort has been di- 
rected to assess potential performance of quasi- ID Field Effect 
Transistors (FETs) based on Carbon Nanotubes 0], 0, 
(CNTs), Silicon NanoWires g) (SNWs), Graphene nanorib- 
bons (GNRs) versus the International Technology Roadmap 
for Semiconductors [5] (ITRS) requirements, both from an ex- 
perimental and a theoretical point of view. However, attention 
has been focused on electrical quantities like I on /I s, sub- 
threshold slope, mobility, transconductance |6), Q, 0, while 
an accurate investigation of electrical noise has been often 
neglected. Although the 1/f noise represents the major noise 
source affecting CNT-FETs performance |9|, [10|, the intrinsic 
shot noise is not only critical from an analog and digital 
design point of view, but can also provide relevant information 
regarding interactions among carriers ITD l[T2l |fl~3), electron 
energy distribution [14] [15] and electron kinetics fl6l . 

Due to the limited device dimensions, even in strong 
inversion only few electrons take part to transport, so that 
drain current fluctuations can heavily affect device electrical 
behavior. Pauli and Coulomb interactions play an important 
role in noise analysis, through fluctuations of the occupation 
number of injected states and fluctuations of the potential 
barrier along the channel. 

From a numerical point of view, a self-consistent solution of 
the electrostatics and transport equations is mandatory in order 
to properly consider such effects. An analysis of this kind has 
been performed for example in double gate MOSFETs [17] 



and in nanoscale ballistic MOSFETs lfl8l . where a strong shot 
noise suppression, mostly due to Pauli exclusion principle, has 
been observed. 

A different approach, based on quantum trajectories within 
the De Broglie-Bohm framework has been instead presented 
in |fl9l , where resonant tunneling diodes have been studied 
and heavy approximations have been adopted in order to 
easily consider electron-electron correlation in the many-body 
problem. 

Actually, a complete understanding of the mechanism of 
suppression of shot noise in CNT and SNW-FETs is still a 
debated issue. Indeed, the significant suppression of current 
fluctuations by more than a factor 100 obtained at low temper- 
ature for suspended ropes 0.4 /im long of single wall carbon 
nanotubes [20 1 has not been supported by a comprehensive 
theoretical analysis. Recent experiments of shot noise in CNT- 
based Fabry Perot interferometers lETI show that by only 
including Pauli exclusion one is able to explain most of the 
dependence of shot noise on the backgate bias, but in some 
bias conditions additional mechanisms of electron-electron 
interaction might be needed to explain the observed noise 
suppression. Theoretical efforts have been mainly addressed 
to model the electrical noise in SNW-FETs, where, within a 
scattering approach with the limitation of excluding space- 
charge effects on electron transmission, Pauli exclusion re- 
duces electrical noise in strong inversion down to 0.6% of 
the full value for a gate overdrive of 0.3 V l22l . whereas an 
interesting increase of noise is observed by including electron- 
phonon scattering processes [23|. 

Here, we present a new method to compute the shot 
noise power spectral density in ballistic CNT and SNW-FETs 
based on Monte Carlo (MC) simulations of randomly injected 
electrons from the reservoirs. In order to consider correlations 
between fermions, an analytical formula for the noise power 
spectral density has been computed by means of a statistical 
approach within the second quantization formalism. The de- 
rived formula has then been implemented in the self-consistent 
solution of the 3D Poisson and Schrodinger equations, within 
the NEGF formalism. 

II. Theory 

The average current in a mesoscopic conductor can be 
expressed by means of Landauer's formula: 



{I) = ihJ ^{M^)] Us(E)-f D (E)}} 



(1) 
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where t is the transmission amplitude matrix for states emitted 
from the source (S) and collected at the drain (D) and fg and 
jo are the Fermi-Dirac statistics of the S and D, respectively. 

The zero-frequency noise power spectral density for a two- 
terminal conductor - the so-called Landauer-Biittiker noise 
formula - reads (24) (25): 

JdE{[j s {l-js) + jnil-jo)} Tv[thth] 

+ [fs(l-jD)+jn(l-js)](Tr[th]-Tr[thth])} , 

(2) 

where is the conjugate transpose of the matrix t. However, 
eq. (|2| holds only if one assumes that fluctuations of the 
potential profile do not occur, i.e. that Coulomb interaction be- 
tween carriers is completely neglected. Actually, the potential 
barrier along the channel fluctuates in time, since randomly 
injected electrons modify the height of the barrier through 
long-range Coulomb interaction, which in turn affects carriers 
transmission and eventually leads to the suppression of the 
drain current fluctuations. 

In order to compute the expression of the power spectral 
density in the general case, we take advantage of the second 
quantization formalism. In particular, at zero magnetic field, 
the time-dependent current operator at the source can be 
expressed as the difference between the occupation numbers 
of carriers moving inward and outward the source contact in 
each quantum channel (24l (ng m and rig m , respectively): 



1(f) 



^ J2 JdE [n + Sm {E,t)-n Sm {E,t)} , (3) 



where 



(E,t) = / d(hu)a+ m (E)a Sm (E + hu)e- 



(E,t) = / d(tkj)b+ m {E)bs m (E + tiuj), 



The operators a Sm (E) and as m (E) create and annihilate, 
respectively, incident electrons in the source lead with total 
energy E in the transverse channel to. In the same way, 
the creation b Sm (E) and annihilation bs m (E) operators refer 
to electrons in the source lead for outgoing states. For the 
CNT case, the channel index rn runs over all the transverse 
modes and different spin, whereas for SNW, it also runs along 
the six minima of the conduction band in the k space. In 
addition, the operators as and b$ are related through the 
unitary transformation 



(E) 



(E)a an (E), 



(5) 



where the scattering matrix s has dimensions (Ns + No) x 
(TVs + No) and Ns and No are the number of quantum 
channels in the source and drain contacts, respectively. In the 
following, time dependence will be neglected, since we are 
interested to the zero frequency case. 

If a) is a many-particle (antisymmetrical) state, the oc- 
cupation number a am (E) in the reservoir a (a — S, D) in 
the channel m can be either or 1, and can be expressed 



as VamiE) = (a^ arn (E)a arn (E)) a . Since we are interested to 
current fluctuations, we need to consider an ensemble of many 
electrons states {| eri), | 0-2)1 1 cr 3 ),••-, ctjv}} and to compute 
statistical averages () s . By assuming no correlations between 
states at different energy or injected from different reservoirs, 
the statistical average of <J am {E) reads 



(cr am (E)) s = {{al n (E)a am (E)) a ) s = j a (E) 



(6) 



In the following, we identify (( ) a ) s with ( ). By means of 
(|5j, we obtain the mean current: 



tf> - 55/* {£<[**<*>]. 



(7) 
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where t' is the drain-to-source transmission amplitude ma- 
trix [ |26) . Since a 2 am — a am , Vm G a and exploiting the 
unitarity of the scattering matrix, the mean squared current 
fluctuation for unit of energy can be expressed as: 

var (I) / e\ 2 , 



AE 
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« r~t 
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(8) 



(4) where [t] is defined as 



PL 



[ttt] 
rt'tfi 



if a = S 
if a = D , 



and r is the reflection amplitude matrix |26|. AE is our energy 
step of choice, i.e. the minimum energy separation between 
injected states. 

Eq. ^ is expressed as the sum of four terms: the first, 
the second and the third terms correspond to the partition 
noise contribution. In particular, the first term is strictly related 
on the quantum uncertainty of the transmission process and 
disappears in the classical limit; the second term is associated 
to the correlation between transmitted states coming from 
the same reservoir; the third term to the correlation between 
transmitted and reflected states in the source lead; the minus 
sign in the second and third terms is due to exchange pairings, 
because of the fermionic nature of the electrons. In particular, 
the second and the third terms provide physical insights on 
exchange interference effects (27l . Finally, the last term rep- 
resents the injection noise obtained as the variance computed 
on the ensemble of current samples. 

According to the Milatz Theorem (28) , the noise power 
spectral density in the zero frequency limit can be computed 
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as 5(0) = lirrif-yo S(f) = lim u ^Q [2/v ■ var(I)), where v is density becomes: 
the injection rate, which can be expressed as: 



v = AE/(2irh). 



(9) 



Eventually, the power spectral density of shot noise at zero 
frequency can be expressed as: 



2 unr(I) 
5(0) = lim - var(I) = lim 4irh ^ (10) 



It is worth noticing that eqs. ( fTOb and ((HJ are not equivalent 
to the Landauer-Buttiker's formula (O, since in eq. ([H} the 
transmission (t, t') and reflection (r) matrix are expressed 
as functionals of the statistics of the occupation of injected 
states from both contacts. In this way we are able to consider 
the fluctuation in time of the conduction and valence band 
edge profiles produced by the random injection through long- 
range Coulomb repulsion, providing a further source of noise 
suppression not included in eq. @. 

Indeed, from an analitical point of view, eqs. (0 and (TTOb 
reduce to eq. (fJJ when transmission and reflection do not 
depend, through Coulomb interaction, on random occupa- 
tion numbers of injected states: in that case we can take 
the terms related to transmission and reflection out of the 
statistical averages in (0. By means of d6) and exploiting 
(a al {E)a 0n {E')} s = f a (E)fp(E') + 6(E-E')6 afj 6 ln [f a (E) 
—f a (E)fp(E')], the fourth term in $8$ becomes: 



(!) dE^[th(B)] 2 nnfs (E)[l-ME)} 

neS 

dEj2 [ttt(E)]J fc / D (^0 [l-fn{E)] (11) 



h 
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since at zero magnetic field t'+t' = tt+. The terms S(E — 
E'), 6 a p and <5/„ are the Rronecker delta. Taking advantage 
of E k enZ P es [^] kp [^t'] pk = Tr [th] Tr [t+ttt t ], 
5(0) becomes: 



<?((,)= lim 4nh V -^l 
y ' AB^O AE 

=^nJ dE i \fs0--fs) + foO-fo)] Tr [t+tttt] 

+ [fs(l-fD) + f D (l-fs)}(Tr [ttt]-Tr [t+tt+t])} 



(12) 



which is Landauer-Buttiker's formula (|2). 

Let us now point out that eq. ( TTOb would also simplify when 
identical and independent injected modes from the reservoirs 
are considered. In this case, t, t' and r are all diagonal, so that 
the second term in ^ becomes negligible. By exploiting the 
reversal time symmetry (t' = t*, where t* is the transpose of 
t) and the unitarity of the scattering matrix, the power spectral 



[ J a=S,D l£a 

- 2 / dE^2([t\ S]U {l-[t\ s . u )a Dl asi) s 
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III. Simulation Methodology 



In order to properly include the effect of Coulomb inter- 
action, we self-consistently solve the 3D Poisson equation 
imposing Dirichlet boundary conditions in correspondence of 
the metal gates, and null Neumann boundary conditions on 
the ungated surfaces which define the 3D domain. Within a 
self-consistent scheme, the 3D Poisson equation is coupled 
with the Schrodinger equation with open boundary condi- 
tions, within the Non-Equilibrium Green's function (NEGF) 
formalism which has been implemented in our in-house open 
source simulator NanoTCAD ViDES 11291 . In particular the 3D 
Poisson equation reads 



V(eV*(f)) = -G>(0+P/fa , 



(14) 



where <f> is the electrostatic potential, pfi X is the fixed charge 
which accounts ionized impurities in the doped regions, while 
p is the charge density per unit volume 



p(f) 



r+oo 

(E) 

J Ei c r-> ^ r- ... 



a=S,D nea 



+ e 



rE t 

7 ^E Y,DOS an (?,E)[l-a an (E)],(15) 

a=S,Dnea 



where Ei is the mid-gap potential, DOS an {r, E) is the local 
density of states associated to channel n injected from contact 
a and r is the 3D spatial coordinate. 

From a numerical point of view, in order to model the 
stochastic injection of electrons from the contacts, a statis- 
tical simulation on an ensemble of random configurations of 
injected electron states from both contacts has been performed. 
In particular, we have uniformly discretized with step AE the 
whole energy range of integration [equations <[8J an d <TT3T>1 - 
Each random injection configuration has been obtained by 
extracting a random number r uniformly distributed between 
and 1 for each state represented by energy E, reservoir a, 
and quantum channel n. The state is occupied if r is smaller 
than the Fermi Dirac factor, i.e. <rs n (E) [<7Dn(E)] is 1 for 
r < fs(E) [f D (E)], and otherwise. 

Self-consistent simulations for a given actual random statis- 
tics in the source and drain contacts have been then performed, 
and, at convergence, the transmission (t, t') and reflection 
(r) matrix have been computed, obtaining an element of the 
ensemble. In particular, for an actual electron distribution in 
the contacts, the Schrodinger equation is solved in order to 
obtain the spatial charge distribution (TOT l along the channel. 
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Then, the latter is included in eq. (TBI and the electrostatic 
potential is then computed and, once convergence of the 
NEGF-Poisson iteration scheme is achieved, the scattering 
matrix is evaluated, and a new sample to be added to the 
noise ensemble is obtained. Finally, the power spectral density 
5(0) can be extracted by means of eqs. (JHJ and (TTOb . From 
a computational point of view, we have verified that 5(0) 
computed on a record of 500 samples, using the energy step 
AE = 5 xl0~ 4 eV, represents a good tradeoff between 
computational cost and accuracy of results 1 30 1 . 

Let us mention the fact that our approach is based on a 
mean field approximation of the Coulomb interaction, and that 
therefore the exchange term is not included. In the following, 
we will refer to self-consistent (SC) simulations when the 
Poisson-Schrodinger equations are solved considering fs and 
fjj in eq. < fT3T >, while we refer to self-consistent Monte 
Carlo simulations (SC-MC), when statistical simulations with 
random occupations crg n (E) and ao n (E) inserted in eq. (Tl3T > 
are used. SC-MC simulations of randomly injected electrons 
allow to consider both the effect of Pauli and Coulomb 
interaction on noise. As a test, we have verified that if we 
perform MC simulations, keeping the potential profile along 
the channel fixed and exploiting the one obtained by means of 
SC simulation, the noise power spectrum computed in this way 
reduces to the Landauer-Biittiker's limit ©, as already proved 
in an analitical way (eq. (fTSt): we refer to such simulations as 
non-self consistent Monte Carlo simulations (non SC-MC). 

IV. Simulation Results 

A. Considered devices 

The simulated device structures are depicted in Fig. [T] We 
consider a (13,0) CNT embedded in Si02 with oxide thickness 
equal to 1 nm, an undoped channel of 10 nm and n-doped CNT 
extensions 10 nm long, with a molar fraction / = 5 x 10~ 3 . 
The SNWT has an oxide thickness t ox equal to 1 nm and 
the channel length L is 10 nm. The channel is undoped and 
the source and drain extensions (10 nm-long) are doped with 
Nd = 10 20 cm -3 . The device cross section is 4x4 nm 2 . 
A p 2 -orbital tight-binding Hamiltonian has been assumed for 
CNTs lOTl . ||32), whereas an effective mass approximation 
has been considered for SNWTs 11331 . 11341 by means of an 
adiabatic decoupling in a set of two-dimensional equations in 
the transversal plane and in a set of one-dimensional equations 
in the longitudinal direction for each ID subband. 

For both devices, we have developed a quantum fully 
ballistic transport model with semi-infinite extensions at their 
ends. A mode space approach has been adopted, since only 
the lowest subbands take part to transport: we have verified 
that four modes are enough to compute the mean current both 
in the ohmic and saturation region. All calculations have been 
performed at the temperature T= 300 K. 

B. DC Characteristics 

In Fig. |2] the transfer characteristics for different drain- 
to-source biases Vds computed performing SC and SC-MC 
simulations are plotted as a function of the gate overdrive 
Vgs - Vth in the logarithmic scale, both for CNT and SNW 



metal gate 




tox=1 nrrif 



metal gate 



Fig. 1. 3-D structures and transversal cross sections of the simulated CNT 
(top) and SNW-FETs (bottom). 



devices. In particular the threshold voltage Vth for the CNT- 
FET at Vds = 0.05 V and 0.5 V is 0.43 V, whereas we obtain 
V th = 0.13 V for Vds = 0.05 V and 0.5 V for the SNW-FET. 
As can be noted, SC and SC-MC simulations give practically 
the same results for CNT-FET, except in the subthreshold 
region where an interesting rectifying effect of the statistics 
emerges in the Monte Carlo simulations for a drain-to-source 
bias V DS = 0.5 V. 

Instead, the rectifying effect is larger for SNW-FET, differ- 
ences up to 30 % between the drain current (J) computed by 
means of SC-MC and SC simulations can be also observed in 
the above threshold regime. In particular, for a gate voltage 
Vgs = 0.5 V and a drain-to-source voltage Vds — 0.5 V, 
the drain current (I) holds 2.42 x 10~ 5 A applying eq. (0, 
and 1.89 x 10~ 5 A applying Landauer's formula (Q~|). Current 
in the CNT-FET transfer characteristics increases for negative 
gate voltages due to the interband tunneling. Indeed, the larger 
the negative gate voltage, the higher the number of electrons 
that tunnel from bound states in the valence band to the drain, 
leaving positive charge in the channel, which eventually lowers 
the barrier and increases the off current (35). 

In the inset of Fig. [2]the average number of electrons inside 
the channel of a CNT and SNW-FET for two different biases 
Vds = 0.5 V and 0.05 V is shown. As can be seen, only very 
few electrons contribute to transport at any give instant, which 
requires us to attently evaluate the sensitivity of such devices to 
charge fluctuations: the smaller the drain-to-source voltage, the 
larger the average number of electrons in the channel, since, 
for low Vds, carriers are injected from both contacts. 

C. Noise 

Let us now focus our attention on the Fano factor F, defined 
as the ratio of the computed noise power spectral density 5(0) 
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gate voltage (V) 



Fig. 2. Transfer characteristics computed for Vos~ 0-5 V and 0.05 V, 
obtained by SC-MC and SC simulations, for CNT and SNW-FET. Full dots 
refer to CNT, empty dots to SNW. Inset: average number of electrons in 
CNT-FETs and SNW-FETs channel, evaluated for V DS = 0.5 V and 0.05 V. 

and the full shot noise 2e(J), F = S(0) /(2e (I)). In Fig. [3] 
the Fano factor for both CNT-FETs and SNW-FETs is shown 
for Vds= 0.5 V as a function of drain-to-source current (/). 

Let us discuss separately the effects of Pauli exclusion alone 
and concurrent Pauli and Coulomb interactions. Triangles in 
Fig. [3] refer to F computed by means of non SC-MC sim- 
ulations on 10 4 samples, while diamonds to results obtained 
by means of Landauer-Btittiker's formula, applying eq. I©. 
As expected the two approaches give the same results for 
both structures. Solid lines refer to S(0) computed by means 
of eqs. (O and ( fTQb and SC-MC simulations, i.e. Pauli and 
Coulomb interactions simultaneously taken into account. 

In the sub-threshold regime ((/} < 10~ 9 A), drain current 
noise is very close to the full shot noise, since electron-electron 
correlations are negligible due to the very small amount of 
mobile charge in the channel. 

From the point of view of eq. (gj, for energies larger than 
the top of the barrier, we have fo{E) -C fs(E) <C 1 and 
the integrand in © reduces to Tr [t t t(£')] fs{E). Instead, 
for energies smaller than the high potential profile along the 
channel, IVt(-E)] <C 1 Vn, m G S, so that we can neglect 
Tr [t+tttt] in ©,"with respect to Tr [tH]. Since f D (E) < 
fs(E), the integrand in © still reduces to Tr [t^t(E)] f s (E). 
The Fano factor then becomes 

S(0) ^ ^JdETr[tH(E)]f s (E) _ 

2e(/> ~ 2e^JdETr[tH(E)}f s (E) V 

In the strong inversion regime instead ({/) > 10~ 6 A), 
the noise is strongly suppressed with respect to the full 
shot value. In particular for a SNW-FET, at (I) ss 2.4 x 
1(T 5 A (V GS - Vth ~ 0.4 V), combined Pauli and Coulomb 
interactions suppress shot noise down to 23 % of the full shot 
noise value, with a significant reduction with respect to the 
value predicted without including space charge effects as in 
Ref. E2, while for CNT-FET the Fano factor is equal to 0.27 




current (|o,A) current (|aA) 

Fig. 3. Fano factor as a function of the drain current {/) for a) CNT- and b) 
SNW-FETs for V DS = 0.5 V. Solid line refers to the Fano factor F obtained 
by means of SC-MC simulations, dashed line (diamonds) applying eq. (5) 
and dotted line (triangles) by means of non SC-MC simulations. 

at (I) « 1.4 x 1(T 5 A (Vqs - V th ~ 0.3 V). Indeed, an 
injected electron tends to increase the potential barrier along 
the channel, leading to a reduction of the space charge and 
to a suppression of charge fluctuation. Note that, by only 
considering Pauli exclusion principle, we would overestimate 
shot noise by 180 % for SNWT ((I) « 2.4 x 10" 5 A) and 
by 70 % for CNT-FET ((/) « 1.4 x 10~ 5 A). 

D. Shot noise versus thermal channel noise 

According to the classical approach for the formulation of 
drain current noise, channel noise is tipically described in 
terms of a "modified" thermal noise, as 5(0) = jSt, where 
St = ^KsTgda is the thermal noise power spectrum at zero 
drain-to-source bias Vds, fcs is the Boltzmann constant, 7 is 
a correction parameter and gdo = {9(I)/dVDs)v DS =o ^ s tne 
source-to-drain conductance at zero Vds- 

Although the classical formulation accurately predicts drain 
current noise in long channel MOSFETs, where 7 is equal 
to 1 in the ohmic region and 2/3 in saturation |28l . it 
underestimates noise in short channel devices. In particular, ex- 
perimental evidences ll36ll of an excess noise in short channel 
MOSFET have been explained in terms of the limited number 
of scattering events inside the channel which is uneffective in 
suppressing the non-equilibrium noise component 11371 . or in 
terms of a revised classical formulation by considering short 
channel effects, such as the carrier heating effect above the 
lattice temperature [38]. 

Actually, it can clearly be seen that non equilibrium trans- 
port easily provides 7 > 1 and that the cause of 7 > 1 is 
simply due to the fact that channel noise can be more properly 
interpreted as shot noise. For example, in the particular case of 
ballistic transport considered here, we can plot 7 as S(0)/St 
as a function of the gate voltage in Fig. Ht. As can be seen, 
values of 7 larger than 1 can be easily observed in weak 
and strong inversion. The strange behavior of 7 as a function 
of the gate voltage is simply due to the fact that one uses 
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Fig. 4. (a) noise power spectral density obtained by SC-MC simulations and 
thermal noise spectral density as functions of the gate voltage for a) CNT-FETs 
and b) SNW-FETs: the considered drain-to-source biases (Vos = 0-5 V, 0.05 
V) are shown in brackets; c) ratio between the noise power obtained by SC- 
MC simulations and the thermal noise density as a function of the gate voltage. 
gdO is the conductance evaluated for Vds = V: g^o = ( gv^ I 




# of electrons/L (nm 1 ) # of electrons/L (nm 1 ) 

Fig. 5. Fano factor as a function of the average number of electrons inside 
the channel per unit length for three different (13,0) CNT-FETs: (A) t ox = 1 
nm, L= 6 nm, (B) t ox — 1 nm, L= 10 nm and (C) tox= 2 nm, L— 10 nm. In 
a) only the effect of the Pauli principle is shown [eq. (5J]; in b) the effect of 
both Pauli and Coulomb interactions is considered. The drain-to-source bias 
Vos is 0.5 V. 

an inadequate model (thermal noise) corrected with the 7 
parameter to describe a qualitatively different type of noise, 
i.e. shot noise. 

E. Effect of scaling on noise 

Let us now discuss the effect of scaling on noise, focusing 
our attention on a (13,0) CNT-FET. One would expect that 
an increase of the oxide thickness would reduce the screening 
induced by the metallic gate, so that the Coulomb interaction 
would be expected to produce a larger noise suppression. 
For example, in the limit of a multimode ballistic conductor 
without a gate contact, significantly suppression of about two 



order of magnitude with respect to the full shot value has been 
shown by Bulashenko et al |39l . 

However, Ref. [39] exploits a semiclassical approach as- 
suming a large number of modes and the conservation of 
transversal momentum, i.e. the role of the transversal electric 
field induced by the gate voltage is completely neglected. In 
our case only four modes contribute to transport, while the 
top and bottom gates of the simulated devices partially screen 
the electrostatic repulsion induced by the space charge in the 
channel on each injected electron, so that a smaller noise 
suppression than the one achieved in Ref. ||39ll can be expected. 

The Fano factor as a function of the average number of 
electrons inside the channel for unit length, computed by 
means of SC simulation and applying eq. (0, for three CNTs 
with different oxide thickness t ox and channel length L is 
shown in Fig. [5^: it shows results for CNT with t ox = 1 nm, 
L= 6 nm (A), CNT with t ox = 1 nm, L= 10 nm (B), and 
CNT with t ox = 2 nm, L= 10 nm (C). Fig. shows the 
Fano factor computed by performing SC-MC simulations and 
applying eqs. <[8J an d (TTOb - As can be seen, if the Fano factor 
is plotted as a function of the number of electrons per unit 
length, as in Fig. [5] curves are very close to one another, and 
effects of scaling are predictable. 



V. Conclusion 

We have developed a novel and general approach to study 
shot noise in nanoscale quasi one-dimensional FETs, such as 
CNT-FETs and SNW-FETs. Our first important result is the 
derivation of an analytical formula for the noise power spectral 
density which exploits a statistical approach and the second 
quantization formalism. Our formula extends the validity of 
the Landauer-Buttiker noise formula [eq. ©], to include also 
Coulomb repulsion among electrons. From a quantitative point 
of view, this is very important, since we show that by only 
using Landauer-Buttiker noise formula, one can overestimate 
shot noise by as much as 180%. The second important result 
is the implementation of the method in a computational 
code, based on the 3D self-consistent solution of Poisson and 
Schrodinger equation with the NEGF formalism, and on Monte 
Carlo simulations over a large ensemble of occurrencies, 
with random occupation of electronic states incoming from 
the reservoirs. As a final note, we show that scaling of 
ballistic onedimensional FETs is expected to weakly affect 
drain current fluctuations, even in the degenerate injection 
limit. 



VI. Acknowledgment 

The work was supported in part by the EC Seventh Frame- 
work Program under the Network of Excellence NANOSIL 
(Contract 216171), and by the European Science Founda- 
tion EUROCORES Program Fundamentals of Nanoelectron- 
ics, through funds from CNR and the EC Sixth Framework 
Program, under project DEWINT (ContractERAS-CT-2003- 
980409). 



7 



References 

R. Martel, T. Schmidt, H. R. Shea, T. Hertel, and Ph. Avouris. "Single- 
and multi-wall carbon nanotube field-effect transistors". Appl. Phys. 
Lett, vol. 73, pp. 2447-2449, 1998. 

J. Guo, A. Javey, H. Dai, and M. Lundstrom. "Performance analysis 
and design optimization of near ballistic carbon nanotube field-effect 
transistors". IEDM Tech. Digest, ppp. 703-706, 2004. 
N. Neophytou, S. Ahmed, and G. Klimeck. "Non-equilibrium green's 
function (NEGF) simulation of metallic carbon nanotubes including 
vacancy defects". J. Comput. Electron., vol. 6, pp. 317-320, 2007. 
Y. Cui, Z. Zhong, D. Wang, W. U. Wang, and C. M. Lieber. "High 
performance silicon nanowire field effect transistors". Nano Lett., vol. 
3, pp. 149-152, 2003. 

International Technology Roadmap for Semiconductors, URL 
http://public.itrs.net/. 

T. D. Yuzvinsky, W. Mickelson, S. Aloni, G. E. Begtrup, A. Kis, and 

A. Zettl. "Shrinking a Carbon Nanotube". Nano Lett., vol. 6, pp. 2718- 
2722, 2006. 

V. Perebeinos, J. Tersoff, and P. Avouris. "Mobility in Semiconducting 
Carbon Nanotubes at Finite Carrier Density". Nano Lett., vol. 6, pp. 
205-208, 2006. 

J. Knoch and J. Appenzeller. "Tunneling phenomena in carbon nanotube 
field-effect transistors", phys. stat. sol.(a), vol. 205, pp. 679-694, 2008. 
Y. M. Lin, J. Appenzeller, J. Knoch, Z. Chen, and P. Avouris. "Low- 
Frequency Current Fluctuations in Individual Semiconducting Single- 
Wall Carbon Nanotubes". Nano Lett., vol. 6, pp. 930-936, 2006. 
J. Appenzeller, Y.-M. Lin, J. Knoch, Z. Chen, and P. Avouris. "1/f 
Noise in Carbon Nanotubes Devices-On the Impact of Contacts and 
Device Geometry". IEEE Trans, on Nanotechnology, vol. 6, pp. 368- 
373, 2007. 

G. Iannaccone, M. Macucci, and B. Pellegrini. "Shot noise in resonant 
tunneling structures". Phys. Rev. B, vol. 55, pp. 4539^1550, 1997. 
G. Iannaccone, G. Lombardi, M. Macucci, and B. Pellegrini. "Enhanced 
Shot Noise in Resonant Tunneling: Theory end Experiment". Phys. Rev. 
Lett., vol. 80, pp. 1054-1057, 1998. 

Ya. M. Blanter and M. Biitttiker. "Transition from sub-Poissonian to 
super-Poissonian shot noise in resonant quantum wells". Phys. Rev. B, 
vol. 59, pp. 10217, 1999. 

T. Gramespacher and M. Buttiker. "Local densities, distribution func- 
tions, and wave-function correlations for spatially resolved shot noise at 
nanocontacts". Phys. Rev. B, vol. 60, pp. 2375, 1999. 
O. M. Bulashenko and J. M. Rubl "Shot noise as a tool to probe an 
electron energy distribution". Physica E, vol. 12, pp. 857, 2002. 
R. Landauer. "Condensed-matter physics: The noise is the signal". 
Nature, vol. 392: pp. 58-659, 1998. 

Y. Naveh, A. N. Korotkov, and K. K. Likharev. "Shot-noise suppression 
in multimode ballistic Fermi conductors". Phys. Rev. B, vol. 60, pp. 
R2169-R2172, 1999. 

G. Iannaccone. "Analytical and Numerical Investigation of Noise in 
Nanoscale Ballistic Field Effect Transistors". J. Comput. Electron., vol. 
3, pp. 199-202, 2004. 

X. Oriols, A. Trois, and G. Blouin. "Self-consistent simulation of 
quantum shot noise in nanoscale electron devices". Appl. Phys. Lett., 
vol. 85, pp. 3596-3598, 2004. 

P.-E. Roche, M. Kodak, S. Gueron, A. Kasumov, B. Reulet, and 

H. Bouchiat. "Very low shot noise in carbon nanotubes". Eur. Phys. J. 

B, vol. 28, pp. 217-222, 2002. 

L. G. Herrmann, T. Delattre, P. Morfin, J.-M. Berroir, B. Placais, D. C. 
Glattli, and T. Kontos. "Shot Noise in Fabry-Perot Interferometers Based 
on Carbon Nanotubes". Phys. Rev. Lett, vol. 99, pp. 156804, 2007. 
H. H. Park, S. Jin, Y. J. Park, and H. S. Min. "Quantum simulation 
of noise in silicon nanowire transistors". /. Appl. Phys., vol. 104, pp. 
023708, 2008. 

H. H. Park, S. Jin, Y. J. Park, and H. S. Min. "Quantum simulation of 
noise in silicon nanowire transistors with electron-phonon interactions". 
/. Appl. Phys., vol. 105, pp. 023712, 2009. 

M. Buttiker. "Scattering theory of current and intensity noise correlations 
in conductors and wave guides". Phys. Rev. B, vol. 46, pp. 12485-12507, 
1992. 

Th. Martin and R. Landauer. "Wave-packet approach to noise in 
multichannel mesoscopic systems". Phys. Rev. B, vol. 45, pp. 1742- 
1755, 1992. 

S. Datta. Electronic transport in mesoscopic systems. Cambridge 
University Press, 1995. 



[27] A. Betti, G. Fiori, and G. Iannaccone. "Statistical theory of shot noise 
in quasi- ID Field Effect Transistors in the presence of electron-electron 
interaction", http://arxiv.org/abs/0904.4274 (2009). 

[28] A. van der Ziel. Noise in Solid State Device and Circuits. Wiley, New 
York, pp. 16 and 75-78, 1986. 

[29] Code and Documentation can be found at the url: 
http ://w w w.nanohub. org/tools/vides . 

[30] A. Betti, G. Fiori, and G. Iannaccone. "Shot noise in quasi one- 
dimensional FETs". IEDM Tech. Digest, pp. 185-188, 2008. 

[31] G. Fiori and G. Iannaccone. "Coupled Mode Space Approach for 
the Simulation of Realistic Carbon Nanotube Field-Effect Transistors". 
IEEE Trans, on Nanotechnology, vol. 6, pp. 475^179, 2007. 

[32] J. Guo, S. Datta, M. Lundstrom, and M. P. Anantam. "Towards Multi- 
Scale Modeling of Carbon Nanotube Transistors". Int. J. Multiscale 
Comput. Eng., vol. 2, pp. 257, 2004. 

[33] G. Fiori and G. Iannaccone. "Three-dimensional simulation of one- 
dimensional transport in silicon nanowire transistors". IEEE Trans, on 
Nanotechnology, vol. 6, pp. 524-529, 2007. 

[34] J. Wang, E. Polizzi, and M. Lundstrom. "A three-dimensional quan- 
tum simulation of silicon nanowire transistors with the effective-mass 
approximation". J. Appl. Phys., vol. 96, pp. 2192-2203, 2004. 

[35] G. Fiori, G. Iannaccone, and G. Klimeck. "A three-dimensional 
simulation study of the performance of carbon nanotube field-effect 
transistors with doped reservoirs and realistic geometry". IEEE Trans. 
Electron Devices, vol. 53, pp. 1782-1788, 2006. 

[36] A. A. Abidi. "High-Frequency Noise Measurements on FET's with 
Small Dimensions". IEEE Transactions on electron devices, vol. 33, 
pp. 1801-1805, 1986. 

[37] R. Navid and R. W. Dutton. "The physical phenomena responsible for 
excess noise in short-channel MOS devices". IEDM Tech. Digest, pp. 
75-78, 2002. 

[38] K. Han, H. Shin, and K. Lee. "Analytical Drain Thermal Noise Current 

Model Valid for Deep Submicron MOSFETs". IEEE Transactions on 

electron devices, vol. 51, pp. 261-269, 2004. 
[39] O. M. Bulashenko and J. M. Rubi. "Shot-noise suppression by Fermi 

and Coulomb correlations in ballistic conductors". Phys. Rev. B, vol. 

64, pp. 045307, 2001. 



